Accelerating expansion of the universe may be caused by inhomogeneities 
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We point out that, due to the nonlinearity of the Einstein equations, a homogeneous approxima- 
tion in cosmology leads to the appearance of an additional term in the Friedmann equation. This 
new term is associated with the spatial inhomogeneities of the metric and can be expressed in terms 

_ 2 

of density fluctuations. Although it is not constant, it decays much slower (as t~~s) than the other 
terms (like density) which decrease as t~ 2 . The presence of the new term leads to a correction in 
the scale factor that is proportional to t 2 and may give account of the recently observed accelerating 
expansion of the universe without introducing a cosmological constant. 
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Measurements of the light curves of several hundred 
type la supernovae FIJ2, H Q and other independent 
observations H H @ convincingly demonstrate that 
the expansion of the universe is accelerating. 

This unexpected result stimulated a number of the- 
oretical investigations. Most explanations suggested so 
far seem to belong to one of three categories: assuming 
a nonzero cosmological constant IfA 1 1 fl . assumin g a new 
scalar field ^quintessence" ) [HEiQl El El 111 El 

tha I20L |2 1l. o r assuming new gravitational physics 

In the present paper we follow a fourth approach, 
namely, we attribute the accelerating expansion to a con- 
sequence of the inhomogeneities present in ordinary mat- 
ter, according to the dynamics described by the usual 
Einstein equations. A similar approach has been sug- 
gested in |28j. We assume a matter dominated and (on 
the average) flat universe (zero pressure, fl — 1, k = 0), 
without cosmological constant or any other exotic con- 
stituent. 

It is natural and usual to apply a homogeneous ap- 
proximation for the metric, since universe is indeed ho- 
mogeneous on the large scale. Explicitly, this means that 
instead of the actual space dependent metric gik{t, r) one 
uses its spatial average, i.e. 



9ik(t) = Jim / d 3 r g lk (t,r) 

V^CO V Jy 



(1) 



The actual, not precisely homogeneous metric gik(t, r) 
satisfies Einstein's equations 

1 87tG 
Rik{{9ik}) - ^9ikR{{9ik}) = —^T T ik ■ ( 2 ) 
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Since these equations are nonlinear and since there are 
strong inhomogeneities on smaller scales, the spatially 
averaged homogeneous metric together with the spa- 
tially averaged matter density will not satisfy Eq.©. 
In other words, inhomogeneities (which, unlike a cos- 
mological constant or quintessence, are unquestionably 
present) induce a correction term in the Friedmann equa- 
tion. Such a procedure has been introduced for the case 
of gra vitational radiation in |29| , applied for an ideal fluid 
in |30| and for the case of a scalar field in [3l| . We demon- 
strate this below in the framework of second order per- 
turbation theory. First order corrections are well known 

[USUI 

We write down the spatial average of the 
second order perturbation equations, which (as they are 
linear in the second order correction of the metric and 
the density) enable us to calculate the spatial average of 
the metric and the density. 

The metric is written in the form 



9jk 



a (0) 
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where the upper, bracketed indices refer to the order of 
the perturbation. Henceforth we use comoving coordi- 
nates. Assuming isotropy, we have 



go^(t) = , 
g a p{t) = -R 2 (t)5 a p 



(4) 
(5) 
(6) 



where the bar over quantities means spatial average (cf. 
Eq.Ql). Especially, we have 



9% = -S«e(R 2 ) {i) (t) 



i = 0,1,2 . 

The use of comoving coordinates implies that 



P = P (0) +P (1) + P (2) , 

T 0a = , 
T a @ = 



(7) 
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Zeroth order quantities are those in a flat Friedmann 
model with matter domination, i.e. 
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First order corrections are given by 



with 
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Its solution sounds 
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where A\ is an integration constant. Putting this into the 
Einstein equation we get an ordinary differential equation 
for RW 
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with the solution 
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;H t 
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Here A, B and C Q are functions depending only on the 
spatial coordinates. A and B are arbitrary (determined 
by the initial conditions) while C a has zero divergence. 
Finally, D a p satisfies a wave equation. Note that the 
assumption of isotropy implies C a = and D a p = 0. 
Henceforth we set A = 0, too. 

Second order corrections of the (0,0) component of the 
Einstein equations have the general structure 



«» =A 2 H-lA lB (») + I L|-(vip(5 ff0( 
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The integration constants A\ , A 2 are to be determined 
from the initial conditions for the second order correction 
of the metric and its time derivative. It is interesting to 
compare Eq. (1221) with the perturbation series of a homo- 
geneous open universe, by considering the curvature term 
as a perturbation in a flat universe: 



(terms linear in g^) 
+ (terms quadratic in gj£ ) = 87rG p*- 2 -* . (17) 

Calculating the spatial average of this equation we have 
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Similarly, for the spatial average of the second order 
corrections of the component of the divergence equation 
we obtain 
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In order to check Ea. (|22|l . we determined numerically 
the second time derivative of the spatial metric from the 
Einstein equations by using an effective pseudospectral 
scheme. The relative error of the coefficients (mainly 
due to third order corrections) is displayed in FigQ] 

For sub-horizon-sized perturbations in the late uni- 
verse the last term of Ea. (|22|l is dominant. The complete 
scale factor up to second order then reads 



R(t) = 
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where 8p = p 
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FIG. 1: Relative error of the coefficients of Eg. 1221 versus 
V (Sp/p) 2 . Solid line: relative error of the coefficient of the t 2 
term, dashed line: relative error of the coefficient of the t 4 ' 3 
term. 

The present value of the fractional density fluctuation 
can be estimated as (provided that linear perturbation 
theory is applicable, for — > 1 fluctuations grow faster) 
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(see e.g. |35|). Here and dec stand for present value and 
value at decoupling, respectively. Thus, the final results 
are expressed in terms of the (large) relative density fluc- 
tuations. Ea. H24|) demonstrates that taking into account 
of density fluctuations is essential in the late universe. A 
perturbative treatment may not even be sufficient. It is 
also seen that the additional term due to inhomogencitics 
is a quadratic function of the time, thus for a long time 
after the Big Bang, accelerating expansion occurs. This 
can be quantified in terms of the deceleration parameter 
which reads 
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Thus, according to second order perturbation the- 
ory the deceleration parameter becomes negative (ac- 
celeration) if Sp/p > y^3/7 w 0.65. Assuming that 



(Sp/p) ~ 1, this happens in our approximation at 
z = sfTJl- 1 « 0.53. 

It is instructive to calculate the effective pressure as 
well. In a homogeneous universe with the scale factor 
iPH we obtain 
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This negative pressure is related to the gravitational at- 
traction of the underlying inhomogeneities rather than 
a cosmological constant. Expressing it in terms of the 
density we get up to second order accuracy the equation 
of state 
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In order to have a tentative comparison with super- 
nova measurements, we calculate the luminosity dis- 
tance cLl versus the redshift z (Hubble diagram) by 
approximating the inhomogeneous universe with a ho- 
mogeneous one, which expands according to the cor- 
rected scale factor. Differences in distance modulus 
m — M = 51og 10 ((ii/10pc) (compared to an empty uni- 
verse, <1l (empty) = z + z 2 /2) are displayed in Fig[3 This 
is to be compared with Figs. 10-11. in |l| or Figs. 4-5. in 
[2j. Note that in the applied approximation curves be- 
longing to higher values of relative density fluctuations 
converge to an upper bound which is roughly identical 
with the uppermost curve displayed. 
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FIG. 2: Differences in distance modulus versus the redshift. 
Different curves correspond to different strengths of relative 



(28) 

density fluctuations y (Sp/p) 2 . 



Finally, we calculate the age of the universe according 
to Eq.J23J. By solving the equation R(t )/R(t ) = H a 
we get 
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tn = 
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Assuming (5p/p) « 1 it yields to = f 3^7 

Some remarks are in order. The averaging procedure 
we used is not quite unambiguous. We might have in- 
cluded e.g. the square root of the spatial metric. In this 
special case we obtain a further positive contribution to 
the coefficient of the t 2 term in the scale factor. Never- 
theless, it is advisable to consider directly the physical 
distance between two points as a more significant quan- 
tity. For points that are not very far from each other, 
one has I = y / — g a/ sAx a Axp, or, inserting the previous 
expansions for g a g , 
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- Ax a Ax = R^y/AxJ. By fixing the 

coordinate distance d a = \J Ax 2 a and averaging over both 
points one obtains for sub-horizon-sized perturbations 



i = dJ m + 
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where the second term indicates accelerating expansion. 

As for the Hubble diagram, a more precise calculation 
would take into account the inhomogeneities of the met- 
ric along light trajectories, too, not only in the global 
scale factor. Such a procedure would need the complete 
second order perturbation term, including the spatial de- 
pendence. Note that a nonperturbative approach to this 
problem has been published in Ref.js3|. As the spatial 
average of the second order term has proved to be large, 
it may happen (if the variance is also large) that the ac- 
tual Hubble diagram broadens to a strip. This would 
mean that the rather scattered measurement points fol- 
low a systematics and are not due simply to experimental 
error. 

In conclusion, we have demonstrated that inhomo- 
geneities essentially influence the time evolution of the 
late universe. Our results may even need further correc- 
tions, since the terms treated perturbatively proved to be 
quite large. Nevertheless, it is remarkable that accelerat- 
ing expansion may result from inhomogeneities, without 
assuming a cosmological constant, quintessence, or mod- 
ified Einstein equations. 
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